Active and Nonlinear Microrheology in Dense Colloidal Suspensions 
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We present a first-principles theory for the active nonlinear microrheology of colloidal model sys- 
tems: for constant external force on a spherical probe particle embedded in a dense host dispersion, 
neglecting hydrodynamic interactions, we derive an exact expression for the friction. Within mode- 
coupling theory (MCT) , we discuss the threshold external force needed to delocalize the probe from 
a host glass, and its relation to strong nonlinear velocity-force curves in a host fluid. Experimental 
microrheology data and simulations, which we performed, are explained with a simplified model. 
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Microrheology is a promising technique providing local 
probes of the dynamics in a complex fluid [l|] . Monitoring 
the motion of a singled-out probe particle embedded in 
a (dense) host fluid or gel, one addresses questions about 
the microscopic origins of the host's complex- fluid behav- 
ior and in particular the link between microscopic mech- 
anisms and macroscopic properties amenable to conven- 
tional rheology. This connection subtly depends on the 
host, probe-bath interactions, and on the applied forces. 
Active microrheology turns this into an advantage, at the 
cost of requiring much better knowledge about the mi- 
croscopic processes [2:]: applying a known forcing to the 
individual particle, one explores the noncquilibrium and 
usually nonlinear regime, providing detailed insight into 
the structure-dynamics relationship, e.g., in cellular en- 
vironments [3] or close to the glass transition 0, S Ejl. 
Laser tweezers, magnetic or surface-chemistry forces [7[ 
provide experimental realizations achieving large forcing. 

The external-force- velocity relations obtained in dense 
suspensions reveal striking nonlinearities, induced by the 
slow relaxation of the host. Leaving the linear-response 
regime, a sudden strong increase in the velocity reveals 
the strength required to pull free the probe from the 
(transient) local neighbor cage, such that force-induced 
motion overrules structural relaxation. Recent theoreti- 
cal progress [1, Q notwithstanding, it remains to under- 
stand the nonlinear friction induced by the slow struc- 
tural rearrangements of host particles. 

Here we develop a theory for active nonlinear microrhe- 
ology in suspensions close to their glass transition, es- 
tablishing the conceptual connection between micro- and 
macro-rheology, when (de) localization of the probe oc- 
curs, and how the structure of the cage is distorted close 
to this yielding point. We start from microscopic equa- 
tions of motion and relate the force-velocity relation of 
the probe, by virtue of an exact Green-Kubo-like formula, 
to a microscopic-force autocorrelation function. This 
can be approximated through non-equilibrium tagged- 
particle density correlation functions, which in turn are 
calculated in the framework of the mode-coupling theory 



of the glass transition (MCT) [l3|. For a hard-sphere 
(HS) suspension with a pulled probe of same size as the 
host particles, we demonstrate that the theory predicts a 
delocalization threshold force that explains the nonlinear 
response seen in experiment and simulation. 

We start from the many-body Smoluchowski equation 
for the noncquilibrium distribution function '^{t) of a sys- 
tem of N Brownian particles (positions r^) and a single 
probe (labeled s), 9t*(t) = Q.^{t). Subjecting only the 
probe particle to a constant, homogeneous force F°^, the 
Smoluchowski operator f2 = fip + Afi reads 

f)= ^ 9,-(fcBm-F0/O-(9.-F"'^)/C., (1) 

i=l,...,N,s 

where AO = —{dg ■ F°^)/(s is the noncquilibrium term 
describing active forcing, and Qq the equilibrium time- 
evolution. We neglect solvent-induced hydrodynamic in- 
teractions and introduce Stokes friction coefficients for 
host (Ci=i,...Ar = Co) and probe (Cs) particles. The Fi^s 
are (potential) interaction forces among the particles. 

To obtain noncquilibrium averages formed with the 
force-dependent Smoluchowski operator, the integration- 
through-transients (ITT) formalism [ll| recasts Eq. ([1]): 

*(i) - *eq - ^ /* dt' exp[m'](F<='' • F,)*,^ , (2) 

assuming equilibrium at i = 0. In particular, the sta- 
tionary friction coefficient ({F"^^), defined via the average 
stationary velocity at given external force 

C(-i;.)t^oo =C(^^.)oo (3) 

is found by using Eq. ([2]) and equating external and in- 
teraction forces on the probe: 

1 

C^Cs + ^ dt (F, exp[f^t(F-)t]F,)eq . (4) 

This formally exact generalized Green-Kubo relation con- 
nects the far-from-equilibrium response to a transient 
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equilibrium-averaged correlation function. ITT achieves 
that all following averages are equilibrium ones, denoted 
by (•) (suppressing the eq subscript). is the adjoint of 
fi, and F'^^ enters non-perturbatively; linear response is 
recovered by neglecting this dependence. 

Following MCT, we assume that force fluctuations are 
governed by collective and probe-particle density fluctu- 
ations, Qq — X]il=i 6xp[*9''i] s-iid 8q = exp[iqrrs]. We 
take that in the thermodynamic limit, the motion of the 
probe has negligible impact on the bulk properties of the 
host, and restrict wave numbers to a discrete grid ne- 
glecting anomalous long distance correlations. Inserting 
a projector P2 oc J2kp QkQp){QkQp o^i both sides of the 
operator exponential in Eq. (|3]), because forces on the 
probe relax by host particle rearrangements and probe 
motion, and splitting four-point density averages into 
dynamical density correlators, 4'k{t) — exp[r2tt]gifc) 
and </>fc(i) = (gl^ exp[ri'''t]p^), we arrive at 
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-r^{t)^^k{t) . (5) 



Sk = (QkQ-k) and — {g%g-k) are the equilibrium 
structure functions describing interactions among probe 
and host particles. 

The probe correlator (j)q{t) is complex- valued, as the 
perturbed operator 17^ is non-Hermitian. This reflects 
that the probe-density distribution is shifted by appli- 
cation of an external force: while in equilibrium it is 
centered around the origin, the average position of the 
probe moves, introducing a complex-valued phase factor 
in 4>q{t). Still, Eq. ([5]) maintains C £ R due to the sym- 
metry ^,(0 = 

Equation ([5]) recasts the problem of calculating the 
probe friction as one of calculating collective and tagged- 
particle density correlation functions. To this end, we 



employ Zwanzig-Mori equations of motion 11 1 



dtcl>t{r) - -u:i„4>Ur) - / dt' mUt - t')dt«j>Ut') , (6a) 



closed by the MCT approximation generalizing Eq. ^ 
to finite wave vectors. 



E i^K>^pKi^Ut)Mt)- (6b) 



Again, the physical idea in the MCT approximation is 
that the friction kernel iTi^{t) relaxes by both probe 
and host density dynamics. The coupling coefficients 
are V^^^ = {qp)S^, w^^pS*^, where uj^gp = 

{qk-gT — iF'^^) ■ p/(s- An analogous set of equations 
holds for (j)q{t). Since the external force acts on the 
probe only, the 4>q{t) are in fact determined by the un- 
perturbed Smoluchowski operator, flo, resulting in the 
standard MCT scenario of glassy dynamics 1^, 12|. This 
describes arrest driven by wave vectors connected with a 
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FIG. 1: Threshold force F™(v) needed to delocalize a hard- 
sphere probe particle in a glass of equally large hard spheres 
with packing fraction (p above the glass transition (dotted 
line) , calculated from MCT within the Percus-Yevick approx- 
imation. Blue circles mark i^™ values used in Fig. [2] Inset: 
corresponding schematic-model result (see text). 



typical host particle radius a. Thus the dimensionless 
parameter measuring the effect of the external force is 
aF°^/{k-QT), the work required to pull the probe over 
that distance in relation to thermal energy. 

The macroscopic counterpart to the friction C, is the 
dispersion viscosity rj measured in bulk flow. Within 
ITT, the analog to Eq. g]) holds for the latter [13]. MCT 
expresses this as a functional only of the host correla- 
tors (j)q{t), while in Eq. ([5|), the probe correlators 4'q{t) 
enter. In linear response close to the glass transition, 
identical scaling laws for both closely link micro- and 



macro-rheology [ij]. For large external forces, this cor- 
respondence breaks: Equations ([6]) for the probe corre- 
lator contain a novel delocalization transition that is ab- 
sent in 4>q{t). A probe arrested in a glassy host suspen- 
sion remains localized in its (deformed) nearest-neighbor 
cage (described by — (j)g{t — > oo) > 0), yielding 
zero average velocity (infinite friction) only below a finite 
threshold F°^. At larger force, the probe is pulled free 
(/q = 0) and attains a steady velocity (finite friction) at 
long times. In the liquid, cages are transient, and a rem- 
nant of the threshold survives as a sudden sharp "force 
thinning" in ({F"''). 

The details of the delocalization transition depend on 
the host properties, which we model now as hard spheres 
using the known numerical MCT results for the collec- 
tive density correlators (j)q{t) within the Percus-Yevick 
S'q-approximation [l^. This model yields a glass transi- 
tion at packing fraction (pc ~ 0.516, where (p = {'iiT/3)ga^ 
with number density g is the only parameter. Fig- 
ure [1] shows our results for the delocalization thresh- 
old force, F°'^, for a probe equal to the host particles 
{Sg = 1 + This threshold F^^'iipc) > is finite at 
ifc, and increases further with increasing density. Note 
that F^:^ = 0(50 /cB^/a), much larger than one might 
intuitively expect. This refiects the strong caging force 
exerted by the set of nearest neighbors that must be over- 
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FIG. 2: Contour plot of the probability distribution f"{r) for 
a localized hard-sphere probe of radius a in a hard-sphere sys- 
tem with same radius a,t ip = 0.52, for external forces acting 
to the right with indicated magnitude in units of ksT/a. 



'"\ 

: (a) 


"1 ' 1 

T ▼ 

1 A 


062 
(p = 0.57 : 




M ' 
a " 


q) = 0.55 
If = 0.50 


♦ 


\ L ° 


If = OAS 


: ■' 










. 't. 

1 Y 








Tl.l 


nl 





[kT/a] 




come before the probe can be delocalized. 

The inverse Fourier transform of is the t oo prob- 
ability distribution for the position of a probe starting at 
the origin; Figure [5] shows our results at packing fraction 
ip — 0.52, slightly above the glass transition, for several 
forces below F^^. F'^^ is taken to be in (positive) hor- 
izontal direction, rendering /*(r) rotational-symmetric 
around this axis. For zero force, the distribution is 
spherical-symmetric and centered around the origin; it 
decays on a length scale of 0.2a, the typical localization 
length for solids dominated by hard-core repulsion. Small 
applied forces mainly shift the center of the distribution 
to a position xg « 0.2a, i.e., they push the probe to the 
"cage wall" without essentially distorting the cage. Close 
to the delocalization threshold, however, /^(r) develops 
a deformed tail extending into the force direction, reduc- 
ing the spherical symmetry to a merely rotational one. 
Interestingly, the tail does not extend along the symme- 
try axis; rather, a "dip" is seen in direction of the applied 
force. For F'^^ > F°^, and f^{r) vanish, indicating a 
delocalized state. 

To discuss probe friction or similar dynamical quanti- 
ties, we need to solve the time-dependent and spatially 
inhomogeneous Eqs. A first step towards this is 

to solve a simplified, "schematic" MCT model. Take 
(^l + J^ (l)'{t)(t){t) dt, with 

dt(l)%t)+uj'4>'{t)+ f m'(t~t')dt'(t)'it')dt' ^0, (7a) 

m%t)^Vs^'*m{t), (7b) 

u!^ = 1 — iF°^. The host correlator 0(i) is set by the 
"Fi2 model" often used to describe glassy dynamics in 
equilibrium [l6| . governed by a separation parame- 
ter e such that e < in the liquid, and e > in the 
glass; e thus measures the host interactions. Vs describes 
the strength of probe-host coupling. Its F'^'^ dependence 
is ignored, neglecting the subtle interplay arising from 
static correlations involving more than one wave vec- 
tor; still, the model recovers the qualitative behavior of 
the force threshold (see inset of Fig. [T]). We expect the 
schematic model to be reasonable for small enough ve- 
locities and/or forces, where the friction is dominated by 



FIG. 3: (a) Probe friction ^ as a function of external force for 
packing fractions as indicated, from simulations of a quasi- 
hard-sphere system (symbols), from Brownian dynamics for 
monodisperse HS (Ref. [l3], open symbols), and from the 
schematic model (solid lines; see text) with ksT/a = 0.05586, 
1/Co = 3.44, and parameters {e,Vs) = (0,6), (-0.045,3.5), 
(-0.058,3), (-0.18,2), (-0.5,1.5), (-0.8,1), top to bottom, 
(b) Open symbols: Experimental force-velocity relations for 
a colloidal suspension, from Ref. 0]. Lines: schematic 
model with (£,«») = (-0.006,46), (-0.008,30), (-0.01,15), 
(—0.05, 6), and vi = V2 = 0.1 — «s/5 (bottom to top), shifted 
along X and y by 0.058 and 296. 



universal features of the transition at F°^. 

To test the simplified model, we performed simulations 
of a slightly polydisperse quasi-hard sphere system un- 
dergoing strongly damped Newtonian dynamics, which 
shows a glass transition at ipc ~ 0.595 [l^. Particles 
(mass m = 1, fceT = 1, radii distributed uniformly in 
[0.9,1.1]) suffer friction with the solvent (Co = 50) and 
random forces obeying the fluctuation-dissipation theo- 
rem. One particle is randomly selected to undergo an 
external force F'^^ until it reaches a distance half the size 
of the simulation box (elongated in the direction of F'^^ 
by a factor of 8). The average probe velocity is measured 
sampling more than 300 independent trajectories, and 
the friction is calculated using Eq. ([3]). All simulations 
were initially equilibrated, except for = 0.62, where the 
system was aged for i^, = 25000. At this density, results 
show little influence of ageing for forces F"'*^ > 35 fceT/a. 

A strong decrease in the dynamical friction C around 
= O{40 kBT/a) seen in the simulation [symbols in 
Fig. EJa)] indicates the force threshold. Fitting e and 
Vs per curve, and two shift factors setting the units, the 
schematic model reproduces this behavior for ip < ipc- 
In the (idealized) glass, it predicts a true delocalization 
transition as C ^ oo for F°^ < F^^, exemplified by a 
e = curve in Fig. EJa). In the simulation, ^ remains 
finite in the accessible window presumably because of 
ergodicity restoring processes ignored here [10| . Employ- 
ing larger Vg, our model also explains recent experiments 
on colloidal systems using larger probes as shown 
in Fig. EJb) . In these velocity- force curves obtained be- 
low the glass transition, the force-threshold signature is a 
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FIG. 4: Probe-particle density correlation function (j}q{t) from 
computer simulation at </9 = 0.55 (left) and from schematic 
MCT (right; fit as in Fig. |3]), real (top) and imaginary (bot- 
tom) parts, q \\ F'^^ corresponds to the position of the main 
peak in the static structure factor S{q). For the simulation, 
aF'^^'/ikBT) = 1, 10, 20, 30, 50, 100, and 250 (right to left). 

steep increase of {vs)oo around « 0.2 pN ^ k-QT/a, 
again reproduced by the model. For too large external 
forces, the model fails, as expected above. 

The virtue of the schematic model is to allow more 
detailed qualitative predictions for the slow nonequilib- 
rium dynamics. This is demonstrated by Fig. 31 where 
we compare the tagged-particle density correlation func- 
tion 0q (i) obtained from the simulation for a wave vector 
q II F'^^ with a magnitude corresponding to the nearest- 
neighbor peak in Sq. The simulation confirms the ex- 
istence of complex-valued correlation functions (for this 
q-direction) as a signature of non-equilibrium, naturally 
arising in our microscopic framework. For q perpendicu- 
lar to the external force, (t)q{t) = (0iq(i))* remains real- 
valued, owing to the rotational symmetry (l)q{t) — (j)'Lq{t), 
and exhibits the two-step decay typical for glass formers 
with an intermediate plateau and a final relaxation sped 
up by the external force. No such clear plateau is seen 
in the figure for q parallel to F°^. For large forces, the 
(t>qit) show pronounced oscillations, quite unexpected for 
a Brownian system, and even stronger in the simulation 
data. 

To summarize, we have developed a microscopically 
founded theory for the nonlinear active microrheology 
close to a glass transition. Starting from the Smolu- 
chowski equation without hydrodynamic interactions, 
and applying approximations in the spirit of the mode- 
coupling theory of the glass transition, we predict the 
probe friction as a function of the external force and of 
the equilibrium host structure. 

The theory predics a finite microrheological force 
threshold needed to delocalize a probe from a glassy host, 



locally melting it. In the dense liquid, this is reflected by 
a strong nonlinear decrease in friction coefficients differ- 
entiating the regimes where cages are either broken by 
slow structural relaxation (for small external force), or 
by large enough applied force. A schematic model cap- 
tures these aspects and allows to fit experimental and 
simulation data quantitatively for not too large external 
forces. 

The force threshold could be related to the existence of 
a yield stress well established for glassy colloidal systems, 
and predicted by MCT for constant-velocity bulk driving. 
It will be promising to study more closely this relation 
and the dynamical behavior of the system close to micro- 
and macro-yielding. 
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